Journal  of  Power  Sources  225  (2013)  269-276 


ELSEVIER 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


P  r*! 

Sbb..«ltS 


Finite  element  thermal  stress  analysis  of  solid  oxide  fuel  cell  cathode 
microstructures 

Sushrut  Vaidya,  Jeong-Ho  Kim* 

Department  of  Civil  and  Environmental  Engineering,  University  of  Connecticut,  261  Glenbrook  Road,  Storrs,  CT  06269,  USA 


HIGHLIGHTS 


►  We  develop  a  framework  for  cathode  microstructures  in  solid  oxide  fuel  cell. 

►  We  develop  novel  heuristic  schemes  to  derive  cathode  microstructures. 

►  Temperature-dependent  material  properties  must  be  considered. 
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Two-dimensional  images  of  solid  oxide  fuel  cell  (SOFC)  cathode  microstructures  (50:50  wt.%  LSM:YSZ) 
are  used  to  generate  three-dimensional  finite  element  (FE)  models.  An  approximate,  heuristic  scheme  is 
developed  to  derive  a  microstructure  of  30:70  wt.%  LSM:YSZ  composition  from  the  original,  real 
microstructures.  The  derived  model  is  validated  by  calculating  three-phase  boundary  (TPB)  and  phase 
surface  area  densities  by  comparing  with  data  available  in  the  literature.  Construction  of  such  derived 
microstructures  will  provide  insights  on  the  effects  of  phase  compositions  on  the  mechanics  of  electrode 
structures.  The  heuristic  scheme  is  not  proposed  as  a  replacement  for  rigorous  approaches  such  as  the 
random  packing  model,  but  rather  as  a  simplified  approach  for  deriving  reasonably  realistic  micro¬ 
structures  of  different  compositions  within  a  limited  range  of  validity.  The  models  are  analyzed  using 
finite  elements  to  estimate  thermal  stresses  and  probability  of  failure  using  Weibull  analysis.  The  effects 
of  temperature-dependent  material  properties  and  phase  volume  fractions  on  probability  of  failure  of  the 
cathode  are  discussed. 

Crown  Copyright  ©  2012  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  significant  research  efforts  focused  on  solid 
oxide  fuel  cells  (SOFCs)  have  been  undertaken  in  the  field  of  power 
sources.  SOFCs  are  capable  of  directly  converting  chemical  energy 
into  electrical  energy  with  high  efficiency  and  low  emissions  [1]. 
Anode,  cathode,  electrolyte,  and  interconnect  wires  are  the  basic 
components  of  SOFCs  [1].  Researchers  have  investigated  many 
different  aspects  of  SOFCs,  such  as  configurations  [2,3],  materials 
[4-8],  component  microstructures  [9-12],  electrochemical 
performance  [11—14],  and  thermal  stress  [14-17].  Although  these 
examples  are  representative  of  the  research  areas  that  have  been 
investigated,  they  are  by  no  means  exhaustive.  Researchers  have 
also  employed  a  variety  of  approaches  to  investigate  problems 
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related  to  SOFCs:  analytical  [3],  experimental  [4-12],  and  compu¬ 
tational  [2,13-17]. 

Research  has  established  that  electrochemical  performance  of 
SOFCs  is  significantly  affected  by  component  microstructure  [11,12]. 
In  contrast,  few  studies  have  investigated  the  effects  of  micro¬ 
structure  on  mechanical  performance  [16].  However,  mechanical 
durability  of  the  SOFC  under  steady-state  and  transient  thermal 
loads  is  a  significant  aspect  of  performance;  it  is  thus  important  to 
understand  the  effect  of  component  microstructure  on  probability 
of  failure  under  thermal  loads. 

In  this  work,  we  perform  finite  element  analysis  (FEA)  of 
thermal  stresses  induced  in  reconstructed  SOFC  cathode  micro¬ 
structures  under  spatially  uniform  temperature  fields.  Three- 
dimensional  finite  element  (FE)  models  of  SOFC  cathode  micro¬ 
structures  are  generated  from  a  stack  of  two-dimensional  micro¬ 
structure  images.  The  models  are  subjected  to  FEA  to  determine 
thermal  stresses  due  to  a  steady-state  temperature  change  from 
room  temperature  up  to  operating  temperature.  The  calculated 
thermal  stresses  are  analyzed  using  the  Weibull  method  to 


0378-7753 /$  -  see  front  matter  Crown  Copyright  ©  2012  Published  by  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.jpowsour.2012.10.054 


270 


S.  Vaidya,  J.-H.  Kim  /  Journal  of  Power  Sources  225  (2013)  269-276 


calculate  the  probability  of  failure.  The  cathode  material  is 
50:50  wt.%  LSM-YSZ.  The  effects  of  temperature-dependent 
material  properties  on  the  probability  of  failure  of  the  cathode 
are  investigated.  A  cathode  model  of  30:70  LSM:YSZ  composition  is 
derived  using  a  heuristic  scheme  to  investigate  the  effect  of  varying 
volume  fractions  on  mechanical  behavior  under  temperature  loads. 

Anode  (50:50  wt.%  NiO:YSZ)  and  cathode  (50:50  wt.%  LSM:YSZ) 
microstructures  have  previously  been  analyzed  for  thermal  stress 
using  finite  elements  by  the  authors  [17].  This  study  extends  the 
work  by  including  detailed  validation  of  the  original  and  derived 
cathode  microstructures,  and  by  considering  the  effects  of  varying 
phase  volume  fractions  on  the  mechanics  of  cathode  microstruc¬ 
tures  under  thermal  stress.  A  simple  heuristic  scheme  has  been 
developed  to  derive  cathode  microstructures  of  varying  composi¬ 
tions  within  a  limited  range  from  the  original,  real  microstructures. 
Construction  of  such  derived  microstructures  will  provide  physical 
insights  on  the  effects  of  phase  compositions  on  the  mechanics  of 
electrode  structures  in  SOFC,  which  in  reality  possess  varied 
compositions  through  the  thickness  direction  to  increase  material 
compatibility  between  electrolyte  and  electrodes  and  to  improve 
overall  cell  performance. 

2.  Image-based  finite  element  microstructure  models 

2.1.  Models  based  on  original  images 

Three-dimensional  (3-D)  FE  microstructure  models  were 
reconstructed  from  41  two-dimensional  (2-D)  cross-sectional 
images  of  cathode  microstructure.  The  2-D  images  were  obtained 
by  Dr.  Scott  Barnett’s  research  group  at  Northwestern  University 
[9,10],  using  focused  ion  beam-scanning  electron  microscopy  (FIB- 
SEM).  Examples  of  the  2-D  images  are  shown  in  Fig.  1.  These  images 
are  of  the  real  cathode  microstructure  having  50:50  LSM:YSZ 
composition. 

In  the  SEM  cathode  images,  white  (pixel  value  =  255)  represents 
LSM,  gray  (pixel  value  =  127)  represents  YSZ,  and  black  (pixel 
value  =  0)  represents  the  pores.  The  original  cathode  image 
dimensions  are  217  pixels  (width)  x  147  pixels  (height).  The  x—y 
(in-plane)  spatial  resolution  between  pixels  is  40.8  nm  and  the  z- 
spacing  between  images  is  53.3  nm.  The  cathode  images  are  thus  of 
size  8.85  pm  (width)  x  6.00  pm  (height). 

Finite  element  modeling  was  carried  out  using  the  commercial 
FE  software  ABAQUS  [18].  This  was  done  by  writing  MATLAB  ® 
[19]  programs  to  generate  ABAQUS  input  files  with  all  relevant 
data  for  running  the  FE  analysis.  The  3-D  FE  model  was  recon¬ 
structed  from  the  2-D  images  using  3-D  finite  elements  (8-node 
brick  elements).  In  order  to  increase  computational  efficiency  for 
3-D  analysis  using  finite  element  method,  we  reduced  the  full 
model  to  a  representative  model  in  which  we  sacrifice  some 
details  of  phase  geometry  but  microstructural  skeleton  that  is 
crucial  to  stress  analysis  remains  almost  unchanged.  To  provide 
validation  data  for  this  present  approach,  we  used  2D  cathode 
microstructure  images  and  calculated  three-phase  boundary  (TPB) 
density  and  phase  surface  area  density  for  the  original,  full-size 


50:50  wt.%  LSM:YSZ  cathode,  as  well  as  the  reduced  50:50 
cathode.  The  data  agreed  reasonably  well  with  those  reported  in 
the  literature  [10,12].  This  quantitative  evidence  may  support  the 
validity  of  the  reduced  model  to  a  limited  range  which  provides 
a  basis  for  mechanical  stress  analysis  that  can  achieve  computa¬ 
tional  efficiency. 

The  original  images  of  the  cathode  were  first  simplified  by 
sampling  pixels  at  regular  intervals  to  reduce  the  image  resolution 
while  retaining  a  sufficiently  detailed  microstructure  for  the 
purposes  of  stress  analysis.  The  simplified  cathode  images  were 
3.51  pm  (width)  x  3.02  pm  (height)  in  size.  The  depth  of  the  voxel 
was  53.3  nm.  The  next  step  in  the  model  generation  process  was 
image  stacking.  A  stack  of  all  the  2-D  ( x—y  plane)  images  was 
created  in  the  z-direction  by  using  a  cell  array  construct  to  arrange 
the  images  consecutively.  An  initially  blank  “buffer”  plane  (of  the 
same  size  as  each  image)  was  then  introduced  between  each  pair  of 
consecutive  images.  This  was  necessary  to  achieve  a  simple  “step” 
variation  in  material  properties  between  corresponding  regions  in 
two  consecutive  images.  The  gaps  between  consecutive  images 
were  filled  by  assigning  one  8-node  brick  element  to  each  voxel. 
Consider  the  nth  and  (n  +  l)th  images,  and  the  buffer  plane 
between  them.  The  voxel  connecting  a  given  pixel  on  the  nth  image 
with  the  corresponding  pixel  location  on  the  buffer  plane  was 
assigned  to  the  element  set  corresponding  to  the  value  of  the  pixel 
on  the  nth  image.  Similarly,  the  voxel  connecting  a  given  pixel  on 
the  (n  +  1  )th  image  with  the  corresponding  pixel  location  on  the 
buffer  plane  was  assigned  to  the  element  set  corresponding  to  the 
value  of  the  pixel  on  the  (n  +  1  )th  image.  Thus,  the  geometry  of  the 
cathode  microstructure  was  reconstructed  in  the  3-D  FE  models, 
with  a  step  variation  in  material  regions  between  consecutive 
images. 

The  volume  fractions  of  phases  in  each  3-D  model  were  calcu¬ 
lated  by  counting  the  number  of  voxels  corresponding  to  each 
phase  (based  on  pixel  value),  and  diving  by  the  total  number  of 
voxels  in  the  model.  Information  concerning  the  material  proper¬ 
ties,  boundary  conditions,  initial  temperature,  temperature  field 
and  required  outputs  (e.g.  principal  stresses)  was  also  specified  in 
the  input  file.  The  input  file  was  then  run  using  ABAQUS.  Various 
free-body  cuts  of  the  50:50  LSM:YSZ  3-D  FE  cathode  model  are 
shown  in  Fig.  2. 

2.2.  Derived  model  and  validation 

In  order  to  study  the  effect  of  composition  on  the  probability  of 
failure  of  the  cathode,  we  derived  from  the  original  images  a  3-D 
model  having  volume  fractions  different  than  those  of  the  orig¬ 
inal  50:50  LSM: YSZ  model.  The  compositions  of  the  models 
included  in  this  study  are  given  in  Table  1  below. 

A  heuristic  scheme  was  developed  to  derive  the  3-D  cathode 
model  of  30:70  LSM:YSZ  composition  from  the  original  50:50 
LSM: YSZ  images.  Boundary  pixel  modification  was  used  to  effect 
volume  fraction  alterations.  This  algorithm  involved  the  identifi¬ 
cation  of  pixels  initially  lying  on  two-  or  three-phase  boundaries  in 
the  original  images,  by  comparison  of  pixel  values  with  nearest  in- 


Fig.  1.  Two-dimensional  SEM  images  of  typical  cross-sections  of  original  50:50  wt.%  LSM:YSZ  cathode  [10]. 
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Fig.  2.  Free-body  cuts  of  the  three-dimensional  reduced-size  FE  cathode  model. 


plane  neighbors.  Boundary  pixels  were  identified  as  black  pixels  at 
corresponding  locations  in  a  new  (initially  empty)  cell  array  (the 
“boundary  array”),  while  interior  pixels  were  identified  as  white 
pixels  in  this  new  boundary  array,  which  was  of  identical  size  as  the 
original  image  stack.  The  locations  of  boundary  pixels  in  the  orig¬ 
inal  image  array  were  picked  out  from  the  boundary  array  by 
identifying  black  pixels  in  it.  The  values  of  pixels  adjacent  to 
boundary  locations  in  the  original  image  array  were  modified 
according  to  heuristic  rules  derived  by  trial-and-error,  until  the 
desired  volume  fractions  of  each  phase  were  obtained  in  the  cor¬ 
responding  3-D  model.  The  boundary  pixel  identification  algorithm 
is  illustrated  in  Fig.  3. 

The  steps  in  the  heuristic  algorithm  are  illustrated  in  Fig.  4, 
which  shows  how  the  algorithm  operates  in  a  step-by-step  manner 
on  cathode  images  from  the  50:50  model  to  derive  the  30:70  model 
by  changing  the  relative  proportions  of  pixels  of  different  phases. 
The  microstructural  characteristics  of  the  derived  30:70  cathode 
model  were  validated  by  calculating  three-phase  boundary  (TPB) 
density  and  phase  surface  area  densities  from  the  stack  of  2-D 
images,  and  comparing  these  values  with  data  available  in  the 
literature  (Tables  2  and  3).  The  TPB  density  and  surface  area  density 
calculation  algorithms  were  first  validated  by  calculating  the  TPB 
and  surface  area  densities  for  the  full-size  50:50  cathode  model,  for 
which  data  are  already  available  in  the  literature  [10,12].  TPB  length 
is  calculated  from  the  2-D  image  stack  using  the  boundary  array 
described  above,  by  identifying  three-phase  boundary  points 
through  pixel  value  comparison  with  nearest  neighbors.  TPB 
density  (pm-2)  is  calculated  by  dividing  the  TPB  length  (pm)  by  the 
total  reconstructed  volume  of  the  corresponding  model  (pm3).  The 
phase  surface  area  is  calculated  from  the  boundary  array  by 
counting  the  total  numbers  of  surface  elements  belonging  to  each 
phase.  Finally,  the  surface  area  density  (pm-1)  of  that  phase  is 
calculated  by  dividing  the  total  surface  area  of  the  phase  (pm2)  by 
the  total  reconstructed  volume  of  the  model  (pm3). 

From  Table  2,  we  can  see  that  there  is  reasonable  agreement 
between  the  TPB  and  surface  area  densities  calculated  by  the 
algorithms  for  the  full  and  reduced  50:50  cathode  models,  and 
those  reported  in  the  literature.  For  the  reduced  30:70  cathode 
model  derived  from  the  original  images,  the  calculated  surface  area 
densities  are  in  closer  agreement  with  the  literature  than  the 
calculated  TPB  density.  This  is  probably  due  to  the  fact  that  the  TPB 
density  calculation  algorithm  makes  use  of  the  2-D  images,  while 


Table  1 

Original  (50:50)  and  derived  (30:70)  LSM:YSZ  compositions  of  cathode  [12]. 


Cathode  model 

Composition 

Volume 

Volume 

Volume 

designation 

(weight 

fraction 

fraction 

fraction 

percentage) 

LSM:YSZ 

of  LSM  (%) 

of  YSZ  (%) 

of  pores  (%) 

Cathode  5050 

50:50 

23.56 

24.48 

51.96 

Cathode  3070 

30:70 

15.81 

37.39 

46.80 

the  values  reported  in  [12]  are  based  on  detailed  analysis  of  the  3-D 
microstructure  data  set  and  on  an  analytical  model  [20].  It  is  worth 
noting  here  that  Vivet  et  al.  [21  ]  have  reported  TPB  densities  as  high 
as  7.2—11.2  pm-2  for  Ni-YSZ  anode  cermet  samples  having  NiO  wt.% 
values  varying  from  45%  to  61%,  with  sample  volumes  varying  from 
349.39  to  957.95  pm3.  Moreover,  Wilson  et  al.  [12]  have  pointed  out 
that  TPB  density  may  be  higher  for  cathodes  than  for  anodes  due  to 
the  finer  cathode  microstructures.  In  this  regard,  the  TPB  density  of 
12.9  pm-2  obtained  for  the  reduced  30:70  cathode  model  may  be 
considered  reasonable,  given  that  the  reconstructed  volume  for  the 
reduced  cathode  is  45.2  pm3  (Table  3). 

3.  Finite  element  thermal  stress  analysis 

3  A.  Study  design 

Finite  element  analyses  of  the  cathode  microstructure  models 
were  carried  out  to  investigate  thermal  stresses  due  to  various 
temperature  fields.  The  effects  of  varying  phase  volume  fractions 
and  temperature-dependent  material  properties  on  thermal 
stresses  and  probability  of  failure  were  investigated.  The  FE  anal¬ 
yses  were  divided  into  different  categories  as  explained  in  Table  4. 
In  each  case,  the  FE  model  was  subjected  to  fixed  boundary 
conditions  (i.e.  all  nodes  on  each  of  the  six  faces  were  allowed 
neither  to  translate  nor  to  rotate).  The  behavior  of  the  model  with 
increasing  temperature  loads  was  investigated  by  subjecting  the 
model  to  eight  different  spatially  uniform  predefined  temperature 
fields  of  magnitude  120  °C,  220  °C,  320  °C,  •••,  820  °C.  In  each 
analysis,  the  initial  temperature  was  specified  as  20  °C  (room 
temperature),  so  that  the  model  was  subjected  to  eight  different 
magnitudes  of  temperature  change  (AT  =  100  °C,  200  °C,  300  °C,  •  •  •, 
800  °C). 

3.2.  Material  properties 

Table  5  lists  the  room  temperature  material  properties  used  for 
YSZ  and  LSM. 

Fig.  5  shows  the  variation  of  the  CTE  of  YSZ  with  temperature 
[22]  and  the  variation  of  the  Young’s  modulus  of  LSM  and  YSZ  with 
temperature  [4].  The  CTE  of  LSM  was  assumed  constant  over  the 
temperature  range  considered.  The  room  temperature  value  of  the 
CTE  of  LSM  (Table  5)  was  used  in  the  FE  analyses. 

3.3.  Stress  analysis 

The  3-D  FE  model  of  the  cathode  (50:50  LSM:YSZ)  is  shown  in 
Fig.  6.  The  model  has  244584  elements  and  395131  nodes.  The  von 
Mises  stress  contour  plots  for  the  50:50  LSM: YSZ  cathode  for 
AT  =  100  °C  and  800  °C  are  shown  considering  temperature- 
dependent  material  properties  in  Fig.  7.  The  stress  values  are  in 
units  of  N  m-2,  i.e.  Pa. 


Table  2 

Full  and  reduced  cathode  models. 
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Model  size  and 
wt.%  LSM: YSZ 

TPB  density 
from  algorithm 
(liirr2) 

TPB  density 
from  literature 
(liirr2) 

LSM  SA 
density  from 
algorithm 
(pm”1) 

LSM  SA  density 
from  literature 
(pm"1) 

YSZ  SA  density 
from  algorithm 
(pm”1) 

YSZ  SA  density 
from  literature 
(pm”1) 

Pore  SA 
density  from 
algorithm 
(pm”1) 

Pore  SA 
density  from 
literature 
(pm”1) 

Full,  50:50 

5.74 

7.35  [10] 

1.83 

1.80  [12] 

3.08 

3.00  [12] 

4.05 

3.90  [12] 

Reduced,  50:50 

8.40 

7.35  [10] 

1.67 

1.80  [12] 

2.57 

3.00  [12] 

3.83 

3.90  [12] 

Reduced,  30:70 

12.9 

6-9  [12] 

1.48 

1.30  [12] 

3.68 

4.39  [12] 

3.79 

4.70  [12] 

SA  =  Surface  area. 


Table  3 

Full  and  reduced  volumes  of  cathode  models. 


Model 

X-Y  pixel 

Z-resolution 

Fleight 

Width 

Number  of 

Volume  of 

resolution  (nm) 

(spacing  between 
images)  (nm) 

(pixels  pm-1) 

(pixels  pm-1) 

images  used 

model  (pm3) 

Full 

40.8 

53.3 

147/6.00 

217/8.85 

242 

682 

Reduced 

40.8 

53.3 

74/3.02 

86/3.51 

41 

45.2 

Fig.  7  shows  that  as  AT  increases  from  100  °C  to  800  °C,  the 
stresses  in  the  cathode  also  increase.  Thermal  stress  is  proportional 
to  the  CTE  and  AT,  and  the  CTE  of  YSZ  increases  with  temperature 
while  the  CTE  of  LSM  is  assumed  constant  over  the  temperature 
range  considered.  Also,  the  stress  contours  show  that  the  stresses 
are  greater  near  the  pores  due  to  stress  concentration.  Similar 
results  were  obtained  for  the  case  with  temperature-independent 
material  properties.  The  cathode  model  with  30:70  LSM:YSZ 
composition  shows  similar  behavior  in  both  cases.  The  effects  of 
temperature-independent  versus  temperature-dependent  material 
properties  on  the  principal  tensile  stresses  induced  in  the  cathode 
are  discussed  in  conjunction  with  probability  of  failure  estimation, 
where  the  effects  of  varying  phase  compositions  on  thermal 
stresses  and  probability  of  failure  are  also  discussed. 

4.  Failure  probability  estimation  using  the  Weibull  method 

Ceramic  materials  exhibit  brittle  behavior  under  tensile  stress. 
Also,  unlike  metals,  they  show  wide  variability  in  tensile  strength 
values  and  follow  a  statistical  strength  distribution.  Thus,  the 
Weibull  method  of  analysis  [2,23],  was  used  to  calculate  the 
probability  of  failure  of  the  SOFC  cathode. 

According  to  the  Weibull  method,  the  survival  probability  of 
a  particular  component  j  under  the  action  of  a  tensile  stress  a  is 
given  by  [2]: 


distribution,  while  the  Weibull  modulus  m  is  the  shape  parameter. 
The  reference  volume  Vo  is  related  to  the  characteristic  strength  oo 
of  the  material.  In  this  study,  the  Weibull  method  was  modified  to 
account  for  the  fact  that  the  cathode  material  is  a  composite  made 
up  of  two  different  components  (LSM  and  YSZ).  The  Weibull 
parameters  used  for  the  ceramic  materials  (LSM  and  YSZ)  are 
shown  in  Table  6.  Room  temperature  values  of  the  Weibull 
parameters  were  used  in  this  study. 

The  results  of  each  stress  analysis  were  post-processed  by 
writing  programs  to  extract  the  three  principal  stress  values 
from  each  ceramic  (LSM/YSZ)  element  in  the  cathode  models. 
These  principal  stresses  were  then  used  to  perform  a  Weibull 
analysis  to  determine  the  probability  of  failure  of  the  cathode  at 
each  AT  value.  Since  the  SOFC  component  materials  are  sub¬ 
jected  to  a  multi-axial  state  of  stress,  the  total  survival  proba¬ 
bility  of  each  ceramic  phase  of  the  cathode  under  the  action  of 
the  three  principal  stresses  (04,  (72,  and  (73)  was  calculated.  The 
principal  stresses  were  assumed  to  act  independently,  and  the 
total  survival  probability  was  calculated  as  the  product  of  the 
survival  probabilities  under  the  action  of  each  individual  prin¬ 
cipal  stress  [2]: 


Also, 


Flere  j  =  cathode,  V,  is  the  volume  of  component  j,  Vo  is  a  char¬ 
acteristic  specimen  volume  (reference  volume)  for  the  material  of 
component  j,  00  is  the  characteristic  strength  of  the  material  of 
component  j,  and  m  is  the  Weibull  modulus  of  the  material.  The 
characteristic  strength  ctq  is  also  the  scale  parameter  for  the 


Table  4 

Finite  element  analyses  of  cathode. 


Case 

LSM 

YSZ 

Temperature 

-dependence 

(LSM) 

Temperature 

-dependence 

(YSZ) 

Temperature 

Linear 

Linear 

None 

None 

-independent 

elastic 

elastic 

Temperature 

Linear 

Linear 

Young’s 

Young’s 

-dependent 

elastic 

elastic 

modulus 

modulus,  CTE 

Here,  j  =  LSM  or  YSZ  for  the  cathode,  and  i  =  1,  2,  and  3.  Only 
tensile  values  of  the  three  principal  stresses  were  used  in  the 
Weibull  analysis.  The  probability  of  failure  of  each  phase  was  then 
calculated  as  follows  [15]: 

Pf  =  1.0  —  Ps(ff,  Vj) 


Table  5 

Room  temperature  material  properties  used  in  FE  analyses  [15]. 


Material 

Young’s  modulus  (GPa) 

Poisson’s  ratio 

cte(io-6oc-1) 

YSZ 

205 

0.30 

10.40 

LSM 

40 

0.25 

11.40 
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.  .-6  Variation  of  CTE  of  YSZ  with  temperature 
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Fig.  5.  (left)  Variation  of  CTE  of  YSZ  with  temperature;  (right)  variation  of  Young’s  modulus  of  LSM  and  YSZ  with  temperature. 


4.1.  Probability  of  failure  analyses 

Since  the  cathode  is  composed  of  two  different  ceramic  mate¬ 
rials  (LSM— YSZ),  the  probability  of  failure  of  the  cathode  was 
calculated  by  extracting  positive  (tensile)  values  of  the  three  prin¬ 
cipal  stresses  from  each  element  in  the  LSM  and  YSZ  element  sets  of 
the  cathode  FE  model,  and  subjecting  these  to  Weibull  analyses. 
This  resulted  in  two  different  probability  of  failure  values  for  the 
LSM  and  YSZ  phases  of  the  cathode,  which  were  combined  into 
a  single  probability  of  failure  value  for  the  cathode  by  assuming  that 
the  cathode  fails  when  either  phase  fails  (or  when  both  phases  fail 
simultaneously).  The  probability  that  both  phases  fail  simulta¬ 
neously  was  calculated  by  assuming  that  the  failures  of  the  two 
phases  are  independent  events,  and  hence  the  probability  of 
simultaneous  failure  of  the  two  phases  is  the  product  of  the  prob¬ 
abilities  of  failure  of  LSM  and  YSZ: 

pcathode  =  pf  (LSM  or  YSZ)  =  Pf(LSMUYSZ) 

=  >  Pfcathode  =  pf(LSM)  +  Pf(YSZ)  -  Pf(LSMnYSZ) 

=  >  Ppthode  =  Pf(LSM)  +  Pf (YSZ)  -  Pf(LSM)Pf(YSZ) 

The  probability  of  failure  (Pf)  values  for  the  LSM  and  YSZ  phases 
of  each  cathode  model  were  calculated  and  combined,  as  described 
above,  at  each  AT  value  (100  °C,  200  °C,  •  •  •,  800  °C)  for  both  the 
cases  described  in  Table  4  (temperature-independent  material 
properties  and  temperature-dependent  material  properties).  These 
values  are  plotted  in  Fig.  8. 

The  Pf  plots  for  both  cathode  models  show  that  the  probability  of 
failure  of  the  cathode  increases  with  increasing  AT  values  (and  hence 
increasing  stresses),  for  both  temperature-independent  and 


Fig.  6.  Three-dimensional  reduced-size  FE  models  of  50:50  (left)  and  30:70  (right) 
wt.%  LSM: YSZ  cathode. 


temperature-dependent  material  properties,  as  expected.  These 
values  may  be  verified  against  the  cathode  Pf  values  reported  by 
Anandakumar  et  al.  [15].  They  investigated  the  probability  of  failure 
of  functionally  graded  SOFCs  subjected  to  spatially  uniform  thermal 
loading.  Although  the  microstructures  of  the  SOFC  electrodes  have 
not  been  considered  in  [15],  the  Pf  values  reported  in  that  work  have 
been  calculated  using  the  Weibull  method,  with  the  same  Weibull 
parameter  values  as  used  in  the  present  study.  Anandakumar  et  al. 
[15]  have  reported  (cf.  Figure  16  in  their  paper)  Pf  values  varying  in 
average  10-12  to  10-5  for  graded  and  layered  LSM  cathodes.  The 
cathode  Pf  values  in  the  present  study  vary  between  4.6  x  10  13  and 
4.4  x  1CT6,  which  is  in  good  agreement  with  those  reported  in  [15]. 

4.1.1.  Effects  of  temperature-dependent  material  properties 

As  seen  from  Fig.  8,  higher  Pf  values  are  obtained  when 
temperature-independent  material  properties  (indicated  by  the  solid 
lines)  are  considered,  for  both  cathode  models.  A  physical  explanation 
for  this  observation  is  suggested  by  the  temperature  variation  of  the 
Young’s  modulus  of  YSZ.  For  YSZ,  E  decreases  from  a  value  of  205  GPa 
at  T=  20  °C  to  a  value  of  147.5  GPa  at  800  °C,  as  shown  in  Fig.  5.  On  the 
other  hand,  when  temperature-independent  material  properties  are 
considered,  the  Young’s  modulus  of  YSZ  has  a  constant  value  of 
205  GPa.  This  large  decrease  in  the  Young’s  modulus  of  YSZ  with 
increasing  temperature  leads  to  lower  stresses  in  the  cathode  with 
temperature-dependent  material  properties  than  with  temperature- 
independent  material  properties.  This  in  turn  leads  to  lower  Pf  values 
with  temperature-dependent  material  properties  than  with 
temperature-independent  material  properties.  This  is  confirmed  by 
the  MPTS  plot  for  the  50:50  cathode  model  shown  in  Fig.  9,  which 
compares  the  maximum  principal  tensile  stress  induced  in  the  YSZ 
and  LSM  phases  of  the  model  for  temperature-independent  and 
temperature-dependent  material  properties. 

Fig.  9  shows  that  the  MPTS  induced  in  the  LSM  phase  for 
temperature-dependent  material  properties  is  lower  than  the  MPTS 
in  the  LSM  phase  for  temperature-independent  material  properties 
over  the  entire  temperature  range.  Similarly,  the  MPTS  induced  in 
the  YSZ  phase  for  temperature-dependent  material  properties  is 
lower  than  the  MPTS  induced  in  the  YSZ  phase  for  temperature- 
independent  material  properties  over  the  entire  temperature 
range.  This  implies  that  the  cathode  Pf  values,  which  are  calculated 
on  the  basis  of  the  tensile  principal  stresses  in  the  LSM  and  YSZ 
phases,  will  be  lower  with  temperature-dependent  material  prop¬ 
erties  than  with  temperature-independent  material  properties. 

4.1.2.  Effects  of  variation  of  phase  volume  fractions 

Fig.  8  shows  the  Pf  values  for  the  two  cathode  models  for 
temperature-independent  and  temperature-dependent  material 
properties.  For  temperature-independent  material  properties,  the 
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Fig.  7.  Von  Mises  stress  contour  plots  (in  Pa)  for  50:50  LSM:YSZ  cathode  considering  temperature-dependent  material  properties:  AT  =  100  °C  (left),  AT  =  500  °C  (middle), 
AT  =  800  °C  (right). 


Table  6 

Weibull  parameters  of  ceramic  materials  (room  temperature  values)  [2]. 


Material 

Weibull 

Characteristic 

Reference 

modulus,  m 

strength,  a0 
(MPa) 

volume,  V0 
(mm3) 

LSM 

7.0 

52.0 

1.21 

YSZ 

7.0 

446.0 

0.35 

Fig.  9.  Maximum  principal  tensile  stress  in  LSM  and  YSZ  phases  of  50:50  LSM:YSZ 
cathode  model. 


models  show  almost  identical  Pf  values  over  the  entire  temperature 
range,  as  indicated  by  the  two  solid  lines  that  coincide  with  one 
another.  This  indicates  that  given  temperature-independent 
material  properties  and  similar  pore  volume  fractions  (51.96%  for 
the  50:50  LSMiYSZ  model  and  46.80%  for  the  30:70  LSM:YSZ 
model),  the  variation  in  the  composition  of  the  model  (in  terms  of 
volume  fractions  of  LSM  and  YSZ)  has  very  little  effect  on  the  Pf 
values.  However,  Fig.  8  shows  that  when  temperature-dependent 
material  properties  are  considered,  the  volume  fraction  of  YSZ 
plays  a  significant  role  in  determining  the  Pf  values,  since  the 
Young’s  modulus  of  YSZ  undergoes  a  large  decrease  with  increasing 
temperature.  Based  on  the  previous  discussion,  we  may  anticipate 
that  lower  stresses  will  be  induced  in  the  30:70  model,  since  it  has 
a  higher  volume  fraction  of  YSZ  than  the  50:50  model.  These  results 
are  indeed  observed  in  Fig.  8.  The  30:70  LSM:YSZ  model  shows 
lower  Pf  values  than  the  50:50  model,  for  temperature-dependent 
material  properties,  in  Fig.  8. 

5.  Concluding  remarks 

This  paper  addresses  finite  element  thermal  stress  analysis  of 
microstructure-based  SOFC  cathode  models.  It  investigates  the 
effects  of  temperature-dependent  material  properties  and  varying 
phase  volume  fractions  of  cathode  materials  on  the  behavior  of  the 
cathode  model  under  temperature  loads.  In  this  paper,  an 
approximate  heuristic  scheme  has  been  used,  based  on  boundary 
pixel  manipulation  of  original  50:50  wt.%  LSM: YSZ  microstructure 
images,  to  derive  a  microstructure  of  30:70  wt.%  LSM:YSZ  compo¬ 
sition.  Construction  of  derived  microstructures  using  such  heuristic 
schemes  may  provide  a  simplified,  approximate  approach  for 
investigating  the  effects  of  phase  compositions  on  the  mechanics  of 
electrode  structures,  and  is  not  intended  to  replace  more  rigorous 
approaches  such  as  the  random  packing  model.  The  critical  findings 
are  described  as  follows.  First,  consideration  of  temperature- 
independent  material  properties  of  the  SOFC  cathode  materials 
results  in  higher  probability  of  failure  values  than  those  obtained 
with  temperature-dependent  material  properties.  Second,  the  two 
cathode  models  studied,  with  similar  pore  volume  fractions, 
showed  the  following  trends.  With  temperature-independent 
material  properties,  limited  variations  in  the  volume  fractions  of 
LSM  and  YSZ  phases  had  negligible  effect  on  failure  probability. 
With  temperature-dependent  material  properties,  the  cathode 
model  with  higher  volume  fraction  of  YSZ  was  found  to  have  lower 
stresses  and  lower  probability  of  failure,  due  to  the  decrease  in  the 
Young’s  modulus  of  YSZ  with  increasing  temperature. 

The  probability  of  failure  of  cathode  Pf  values  in  the  present  study 
vary  between  4.6  x  10~13  and  4.4  x  10-6,  which  is  in  good  agreement 
with  those  reported  in  [15]  varying  in  average  10  12  to  10-5  for 
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graded  and  layered  LSM  cathodes.  The  present  work  did  not  consider 
the  effects  of  residual  stresses,  which  can  be  an  interesting  topic  in 
the  future. 
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